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cr: Abstract 

Models of the mind are based on the idea that neuron microtubules 
can perform computation. From this point of view, information process- 
^ ' ing is the fundamental issue for understanding the brain mechanisms that 

produce consciousness. The cytoskeleton polymers could store and pro- 
cess information through their dynamic coupling mediated by mechanical 
energy. We analyze the problem of information transfer and storage in 
brain microtubules, considering them as a communication channel. We 
discuss the implications of assuming that consciousness is generated by 
, the subneuronal process. 



1 Introduction 

FT 

^ . In recent years many papers have addressed the problem of developing a 

theory of mind [1-13]. R. Penrose and S. Hameroff developed a quantum model 
of the mind considering the cytoskeleton of neuron cells as the principal com- 
ponent that produces states of mind or consciousness [2,3]. In their model the 
microtubules (MTs) perform a kind of quantum computation through the tubu- 
lins. Tubulins are proteins which form the walls of the MTs. They claim that 
the tubulins work like a cellular automata performing that kind of computation. 
In this way, the walls of the MT could be able to store and process informa- 
tion by using combinations of the two possible states (a and f3) of the tubulins. 
The MT interior works as an electromagnetic wave guide, filled with water in 
an organized collective state, transmitting information through the brain. A 
gelatinous state of water in brain cells, which was observed by [13], could boost 
these communication effects. 
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Using a different approach, Tuszynski et al. [6-8] model the biophysical as- 
pects of the MTs considering the following questions: What kind of computing 
do microtubules perform? How does a microtubule store and process informa- 
tion? In order to analyze these questions they use a classical approach, studying 
the basic physical properties of the MTs as interacting electric dipoles. 

According to [6-8,14-17] each tubulin has an electric dipole moment ~p due to 
an asymmetric charge distribution. The microtubule is thus a lattice of oriented 
dipoles that can be in random phase, ferroelectric (parallel-aligned) and an 
intermediate weakly ferroelectric phase like a spin-glass phase. It is natural to 
consider the electric field of each tubulin as the information transport medium. 

Therefore, the tubulin dimers would be considered the information unit in 
the brain and the MT sub-processors of the neuron cells. Therefore, to know 
how MTs process information and allow communication inside the brain is a 
fundamental point to understand the mind functions. 

In this work we derive some results which were not explicitly obtained in 
[6-8] and extend the ideas introduced by [6,16] using the point of view of the 
information theory. We analyze the problem of information transfer and stor- 
age in brain microtubules, considering them as a communication channel. The 
electric field is the mediator of each communicator entity. We discuss the im- 
plications of assuming that the consciousness is generated by the microtubules 
as sub-neuronal processors. 

2 Biophysical Aspects of the Microtubules 

The cytoskeleton has a dynamic structure which reorganizes continually as 
the cells change their shape, divide, and respond to their environment. The 
cytoskeleton is composed of intermediate filaments, actin filaments (or microfil- 
aments), and microtubules. The filaments and the microtubules are mutually 
connected and form a three-dimensional network in the cell. There are many 
papers [6-8] showing that the cytoskeleton is the main component which orga- 
nizes the cell, mediates transport of molecules, organelles, and synaptic vesicles. 
The cytoskeleton possibly receives signals from the cellular environment medi- 
ated by the membrane of proteins and participates in signal transmission to the 
neighborhood of the cell [16,17]. 

Microtubules are hollow cylinders whose exterior surface cross-section di- 
ameter measures 25nm with 13 arrays of protein dimers called tubulins. The 
interior of the cylinder contains ordered water molecules which implies the ex- 
istence of an electric dipole moment and an electric field. The MTs represent a 
dipole due to individual dipolar charges of each tubulin monomer. The micro- 
tubule dipole produces a fast growth at the plus end towards the cell periphery 
and a slow growth at the minus end. The MT polarity is closely connected 
with its functional behavior which can be regulated by phosphorylation and 
dephosphorylation of microtubule-associated protein (MAP) [6-8,14-17]. 

Guanosine triphosphate molecules (GTP) are bound to both tubulins in the 
heterodimer. After polymerization, when the hctcrodimer is attached to the 
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microtubule, the GTP bound to the /3-tubulin is hydrolyzed to the guanosinc 
disphosphate (GDP). On the other hand, the GTP molecule of the a-tubulin 
is not hydrolyzed. The microtubules present a calm dynamic instability which 
are their principal feature [6-8] . [?] 

Many models of conformation (and polarity energy) of the microtubular 
protofilamcnt were developed. These models describe the behavior of the pulses 
generated by the free energy in the GTP hydrolysis. The pulses propagate along 
of the MTs through an elastic coupling or through electric field propagation 
between tubulin dimers [5-8,14,15]. The overall effect of the surrounding dipoles 
on a site n can be modelled by the double-well quartic potential [6-7] 



where u n represents the dimcr conformational change on the n-th protofilament 
axis coupled to the dipolc moment. A and B are parameters of the model, 
where A is dependent of the temperature by A = a(T — T c ), T c is the critical 
temperature, and B is a positive parameter independent of the temperature 
[6,8]. In figure 1 we plot the effective potential in terms of u n . 



Our assumptions lead us to reconsider this model taking into account the 
Information Theory to calculate the storage and transference of information 
along the MT. The information is mediated by the electric field propagating in 
the cellular medium. This propagation of energy can provide a communication 
channel. 

3 Communication Channels 

The Shannon entropy of a random variable X is defined by [18]: 



V{u n ) = -jul + 



(1) 



Fig. 1 - Double well quartic 

potential model with a 
potential barrier \A 2 /2B\. 




(2) 
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where p(xi) is the probability of the outcome Xi. This definition describes the 
amount of physical resources required on average to store the information being 
produced by a source, in such a way that at a later time the information can be 
restored completely. 

If we want to send a message X through a noisy channel, that message 
can be subjected to a loss of information. To correlate a sent message X with 
a received message Y we have to calculate the mutual information I(X : Y) 
between them. The mutual information concept gives us how much knowledge 
we obtain from a message X given that we have received Y. It is defined by 
[18,19] 

(I(X : Y)) = (I(X)) - (I(X\Y)) = (I(Y)) (I(Y\X)) (3) 

and 

(i{x\y)) = - E E p^y*) lo sp(^iyfc), (4) 

i k 

where p(xi\y k ) = p(y k , Xi)/p(y k ). 

Nevertheless, by using a binary code to send a message M, compressed by 
procedure C that minimizes the use of bits in that codification, any receiver 
of M, using a decoding procedure D, must to be able to get all information 
associated to M. 

Consider a symmetric memoryless channel 1 N with a binary input Ai n and 
a binary output A ou t ■ For n uses of the channel, the procedure C encodes the 
input message M such that C" : {1, 2 nR } — * Ai n and D decodes the output 
such that D n : {1, 2 nR } — ► A out , where R is the rate of the code (the number 
of data bits carried per bit transmitted) [19]. Therefore, if X is the encoded 
message M through the procedure C, Y is the received message, and D is the 
decoding procedure for Y, then the probability of error is defined by 

p(C n .D n ) = m&xp(D n (Y) ^ M\X = C n {M)). (5) 

AI 

The principal problem of the information theory is to determine the maxi- 
mum rate R for a reliable communication through a channel. When p(C n .D n ) 
— > for n — > oo, the rate R is said achievable. According to Shannon's theorem, 
given a noisy channel N, its capacity f2(iV) is defined to be the supremum over 
all achievable rates for this channel. That is 

n(N)=max((I(X:Y))), (6) 

P(Xi) 

where the maximum is taken over all input distributions p(xi) of the random 
variable X, for one use of the channel, and Y is the corresponding induced 
random variable at the output of the channel. 

Equation J^J allows us to calculate the transference of information among 
many physical systems. The transfer of energy may include the transfer of elec- 
trostatic energy, energy of low frequency oscillating fields, energy of light, energy 

1 The memoryless channel is the one that acts in the same way every time it is used, and 
different uses are independent of one another. 
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of vibrations, etc. Molecules can contain energy in the chemical bonds, in the 
excited electron states, in the conformation states, etc. A common measure of 
the interaction leading to cooperative behaviour is the information transference. 
The electromagnetic field can transfer information through the environment 
among the systems like a communication channel. 

4 Information Processing in Microtubules 

Many features of the cytoskeleton support the idea that microtubules can 
perform computation and store information. According to [6] the charge separa- 
tion of the MTs is wide enough to store information. Due to its dynamic coupling 
the information can be stored as mechanical energy and chemical events. 

Changes in the opposite direction can be favorable to the SG phase over the 
F-phase. This change could switch from the growth mode to operational behav- 
ior. Our focus is this operational mode. Information processing is addressed by 
[1-5] considering the highly specialized nature of the functional proteins on the 
microtubules. 

4.1 Information Storage in Microtubules 

The tubulins form a dipole moment net and therefore are sensitive to exter- 
nal electric fields. Some papers use physical models such as spin net to describe 
the behavior of the dipole moment net [6-7,20]. According to those models mod- 
els, all tubulins are oriented to the same direction at low temperature (~ 200K) 
and the units of the system are organized (figure 2). In this case the system is 
in ferroeletric phase (F). At high temperatures (~ 400K), the system is in the 
paraelectric phase (P) and the polarity of the tubulins are completely disorga- 
nized (figure 3). However, there is a critical temperature T c in which occurs a 
phase transition between F and P, that is, between order and disorder. At this 
phase transition emerges a new state known as spin-glass phase (SG) (figure 4). 
There are some theoretical models trying to estimate this critical temperature. 
One of them estimates the critical temperature around to 300K which is near 
to the human body temperature [7,8]. 



Fig. 2 - schematic picture for 
F-phase 
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Fig. 3 - schematic picture for 
P-phase 



Fig. 4 - schematic picture for 
SG-phase 

We analyze the propagation of informatin along MTs considering the above 
phases. Assuming an energy approximation dependent on the mean polarity 
described by the Landau theory of phase transitions, the total energy can be 
given by [6,16,21] 

£=(|p 2 + ^P 4 )^o, (7) 

where p represents the continuous variable for the mean polarization at each site, 
and Ng is the total number of sites. The parameter a has a linear dependence 
with the temperature a = a(T - T c ), where 200K < T c < 400 K and b > 
[6-8]. E will be minimized by p = for T > T c and by p = ±y/-a{T - T c )/b 
for T < T c . We use the Boltzmann distribution g(p) to weight the energy 
distribution as a function of the mean polarity 

g(p) = Z- 1 e X p(-(3E), (8) 

where /3 _1 = kT, Z is the normalization, and k is the Boltzmann constant. 
Subtituting (JSJ into Q we get 

Because p is a continuous variable, we need to use the continuous counterpart 
of J2j in order to calculate the information mean value of the system. Replacing 
p(x) by g(p) in @ we obtain the following expression for the information storage 
capacity: 

. . InZ a(T — T c ) . 9 . b , ,„„. 

The average of p over the whole MT, considering all domains, is obtained from 

/oo 
g{p)p n dp. (ii) 
-OO 
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Using 1|10[) we can plot the information capacity against the temperature for 
some parameter values. 



Fig 5. - Storage information capacity of MT when 
a = fe = 0.5. 



Fig. 6 - Storage information capacity of MT when 
a=0.5 and b = 50. 
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Fig. 7 - Storage information capacity of MT when 
a = 0.05 and b = 50. 



These graphs corroborate with the results of [6], which show that, probably 
at physiological temperature, we can have a mode of information storage in MTs. 
This is the most important feature for finding another subunit of information 
processing inside the brain. It could show us new perspectives for cognitive 
aspects. 

However, according to these graphs, the maximum information storage is 
obtained at the spin-glass phase, therefore we need to make some assumptions. 
In this phase, there are domains with many energy levels which can store infor- 
mation. The interaction among the domains due to the electric field generated 
by the oscillating dipoles must be considered. This electric field is emitted to 
the neighbouring area producing many channels among the domains in MT. 

4.2 Microtubules as a Communication Channel 

Given the capacity of information storage of MTs, the issue now is to know 
whether there exist some kind of information processing on them. To study any 
kind of processing, it is necessary to describe how the information is stored in 
the MT walls. That is, we need to understand how the information propagates 
along the MT. We saw that the SG phase has the maximum capacity of infor- 
mation storage. Therefore, we will restrict to this phase in order to describe 
the interaction, or communication, among the domains. Here, we are assuming 
that the electric field generated by the MT dipoles is the main mediator which 
allows the communication among the domains. 

The graphs of the previous section show that near to the critical temperature 
T c the information capacity has the maximum capacity of storage. Following 
[6] , we assume in this phase a partition of lattices by local domains (see figure 
4). Therefore, the previous prescription is valid only on the local domains. In 
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this way, a domain j has a polarization pj with probability gj(pj)- If we make 
these assumptions, the total probability is given by 



9 = Y[gj(Pj), (12) 

where r is the number of domains [6] 

As a consequence of (|12|) we have for the spin-glass phase 

3 

with j in the set of domains. 

Now, we need to calculate the amount of information transferred through 
the channels among the domains [6,16]. The domains will communicate only if 
they interact. If we consider two domains, the communication is mediated by 
the electric field interaction between them. In order, to calculate the capacity 
of this communication channel, we use the mutual information concept. 

From © we know how much information is transferred from an event Xk (of 
an ensemble X) to another event yj (of an ensemble Y). The term p(xk\yj) im- 
poses the dependence among the systems. Assuming a Boltzmann distribution, 
we want to know the dependence between the domain k, with polarization pk, 
and the domain j, with polarization pj. This dependence is described by the 
distribution g(pk\pj) which imposes a connection between the domains. Fol- 
lowing [16] we will assume that the output energy is expressed as a function of 
the electric field energy and of the mean polarization energy. Therefore, in the 
thermodynamic equilibrium, we have the average of the output energy E out of 
a domain j as 

{E° ut ) = (Ef gnal ) = (E^ ow ) + (Ej msy ), (14) 

where E^ 9nal is the energy of the coherent signal, Ej msy is the noisy energy, 
and Ej low is the energy of the flow along the system. The energy Ej W is 
responsible for the interaction between the domains. Therefore, supposing that 
the domain j emitts Ej low , we can express the dependence of a domain k as 

g(p k \ Pj ) = Z' 1 exp [-(3{E k + Ef ow + £™)] , (15) 

where Ek is the correspondent energy of the domain k. 

The information entropy depends on the amount of energy in the system 
and on the noisy energy Ej msy . Therefore, the energy of the noise is given by 
[16] 

{E n msy{Tn)) = Z -^MMp))+B(p% (16) 

where A = a(T n - T c ) /2kT n and B = b/AkT n . 

To evaluate the communication channel capacity, each domain is approxi- 
mated by a unique dipole. The information transference will be mediated by a 
radiation of the electric field in the equatorial region of an oscillating dipole. 
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Using the complex Poynting vector and taking the real part, we get an 
expression for the mean value of the flow of energy Ej low . The amount of energy 
absorbed by the oscillating charged units depends directly on its effective cross- 
section and on the intensity of the flow of energy. We can calculate it considering 
the radiation flow of energy through the cross-section D over a spherical surface 
of radius R, where D is a rectangle whose sides are x and z. Hence the expression 
for the flow of energy towards the dipole axis is given by 



(EP° W ) = 2n 2 Sj arcsin 



2R X 



z 1 / z 
2R~ T i\2R z 



(17) 



where Sj = p^rfev^, such that Vj is the dipole frequency, e is the permitivity, 
rj is the permeability of the medium, R x and R z are the perpendicular distances 
from the dipole to the z and x sides, respectively [16]. 

Using JSJ and (0J, and the previous relations, we can derive an expression 
for the capacity of communication between two domains. The channel between 
two domains j and k will be denoted by Njk, hence 

%) = (W-(^I^))> (18) 

where (I(Ek)) is given by an expression similar to l|10|) . The conditional infor- 
mation (I(Ek\Ej)), for a specific polarity pj, can be calculated by a continuous 
version of Q), that is 

(I(E k \E d )) = ^~PJ 9(Pk, Pi) (E k + Ef anal ) dp k . (19) 

Through those calculations we can infer that there is an inter-dependence 
among the domains in the SG phase. Each domain communicates to other 
domain the value of its polarity. It transforms the MT in a net of communication 
units (in this case the units are the domains - see figure 8). Besides, as each 
domain has a particular polarity, in the context of the information theory, we 
can interpret each polarity representing a type of symbol. It would build a 
kind of alphabet along the whole MT, where each domain represents a letter. 
However, as the polarity p is a continuous variable, the change of a letter in 
another would be also in a continuous way and not in a discrete way 



Fig. 8 - A representation of the 
communication between domain j 
and domain k accomplished by the 
electromagnetic field Ej low along 
the walls of MT. 
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Considering the case x = z, we plot the capacity of information transference 
between each domain as a function of the distance and frequency. 



Fig. 9 - Communication channel capacity: 
frequency Vj x distance R z when T ~ 300K. 



Fig. 10 - Communication channel capacity: 
frequency Vj x distance R z when T ~ 1QQK. 
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Fig. 11 - Communication channel capacity: 
frequency Vj x distance R z when T ~ 600K. 

The graphs show that the best conditions to have a communication among 
the domains are at the physiological temperature, with frequency of the confor- 
mational changes of the tubulin dimer protein around to 10 12 s _1 . The relative 
permitivity and permeability in the neighbourhood of the oscillating units is as- 
sumed to be 1 [14-16]. The distance z between the protein molecules is adopted 
to be around to 1/xm — 0.1/zm. At 300K the information transference is su- 
pressed over a distance R z equal to 0.1/xto, and frequency around to O.lTHz 
(figure 9). At a distance smaller than 0.1/um the communication starts to be- 
come independent of the frequency. Finally, at a distance greater thanO.l^m 
the high frequency of the electric field plays a fundamental role in the transfer 
of information. For the other regimes of temperature, the system is not in the 
SG-phase and the graphs show the loss in performance (figures 10 and 11). 

According to [22] , biological molecules with dipolar vibrational activity could 
manifest a quantum coherent mode. That systems could have some isolating 
effect from thermal environments. The frequency range of that quantum mode, 
(also known as Frohlich systems) is around 10 11 s _1 to 10 12 s _1 [1,4]. Therefore, 
the high frequency regimes obtained above can work not only to perform a 
communication along the MT but also to maintain some quantum coherent 
regime 2 . 

5 Conclusions 

This work confirm the results of [1-8] which consider microtubules as a 
classical subneuronal information processor. Through the information theory 
we calculate the information capacity of the MTs. Utilizing models of [1-8] we 
estimate that the favorable conditions for storage and information processing 

2 Somc papers show that the tubulin vibration frequency is in this regime [16,21], 
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are found at temperatures close to the human body. These results corroborate 
the possibility of communication among the domains (where each energy level 
corresponds to some kind of symbol). This communication is mediated by the 
dipole electric field, and this interaction is necessary to describe some processing 
or computing on MT. Through this communication, each domain (or symbol) 
presents some dependence with another. Therefore there are storage as well as 
processing of information associated to the dimers. Besides, from the informa- 
tion theory point of view, the formation of domains creates some redundancy 
for storage or representation of these symbols. This redundancy is important 
for error correction and information protection. However, some points still need 
further investigations. To mention at least two, (1) the direction of the propaga- 
tion of the information under the influence of the environment is an interesting 
point to be analyzed, (2) according to [1-5] there is some water ordination inside 
MTs which could increase the quantum processes in MTs. These points deserve 
to be analyzed using the information theory point of view. 
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